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Abstract. Since the discovery by Hale in the early 1900s that sunspots harbor strong 
magnetic field, magnetism has become increasingly important in our understanding of 
processes on the Sun and in the Heliosphere. Many current and planned instruments 
are capable of diagnosing magnetic field in the solar atmosphere. Photospheric magne- 
tometry is now well-established. However, many challenges remain. For instance, the 
diagnosis of magnetic field in the chromosphere and corona is difficult, and interpreta- 
tion of measurements is harder still. As a result only very few measurements have been 
made so far, yet it is clear that if we are to understand the outer solar atmosphere we 
must study the magnetic field. I will review the history of solar magnetic field measure- 
ments, describe and discuss the three types of magnetometry, and close with an outlook 
on the future. 
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The importance of magnetic field in astropliysical processes has long been recognised. 
Magnetic field is found in all parts of the universe and on all scales: planets, stars, galax- 
ies, accretion disks, etc. Even in cases where magnetic field was initially deemed too 
weak to significantly affect plasma dynamics, it has often turned out to be an important 
factor in the evolution of the system (e.g., in accretion disks, Balbus & Hawley_1991 ) . 
Nowadays a thriving field of research in astrophysics is centered around the life-cycle of 
magnetic field. How is it created? How does it interact with the plasma it is embedded 
in? How is it destroyed? The diagnosis of magnetic field is now spreading throughout 
all parts of astrophysics, e.g., Zeeman-Doppler imaging of stars ( Semel||1989 1 and ex- 
oplanets interacting with their host star ( Shkolnik et al.|2003 1. Indeed, magnetic field 
has now been inferred to be present or even measured in many astrophysical objects. 

Magnetic fields feature perhaps most prominently in solar physics. The definite 
discovery that sunspots have strong magnetic fields was made by Hale ( |1908 1, who 
used a primitive spectro-polarimeter to observe circular polarization in spectral lines 



resulting from the well-known Zeeman effect (Zeeman 1896). A more detailed pub- 



lication ( |Zeema n 1897 a|) was rep roduced in its entirety that same year in The Astro- 
physical Journal ( )Zeeman|1897b) . These discoveries spurred a century of advancement 
spectroscopy, polarimetry, instrument development, and research of solar magnetism. 
It is worth noting that the first detection of widening of spectral lines in sunspots, an ef- 
fect now unambiguously attributed to magnetism, was first reported by Lockyer ( |1866 1, 
more than 30 years prior to the laboratory discovery of the Zeeman effect, and shown 
by [Young] ( 1883 1 in a graphic illustration of an observation dated back to 1870. This 
subject is covered in detail in the excellent review by |del Toro Iniesta ( 1996 1. 
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Figure 1. The Maunder butterfly diagram showing the latitude variation of active 
regions over the solar cycle. 



On the cutting edge are instruments like the Helioseismic and Magnetic Imager 



(HMI, Schou et al. 2012) on the Solar Dynamics Observatory (SDO, Pesnell et al. 



2012[ l, which provides us with the full-disk vector magnetic field in the solar photo- 



sphere several times per hour, and the Spectro-Polarimeter (SP) of the Solar Optical 
Telescope (SOT, Tsuneta et al.|2008 ) of Hinode ( Kosugi et al.|2007 l, which combines 
unparalleled spatial resolution with high sensitivity to make the most refined measure- 
ments of solar photospheric magnetic field to date. 



2. Diagnostics 

Broadly speaking there are two types of diagnostics available to us: qualitative and 
quantitative. The magnetogram, neither qualitative nor quantitative, is a popular diag- 
nostic that occupies the space in between. It is often used because of the difficulties 
associated with performing quantitative diagnostics. 

2.1. Proxy- magnetometry 

Proxy-magnetometry is the technique of determining the locations of magnetic fields 
through a change in intensity. Sunspot maps represent the first application of proxy- 
magnetometry. Over half a century prior to the discovery of the magnetic nature of 
sunspots, the 1 1-year solar cycle was discovered by |Schwabe ( |1844| ), who noticed peri- 



odic behavior after 17 years of observations. Sunspots also emerge closer to the equator 
later in the cycle. Known as Sporer's Law, this feature of the solar magnetic dynamo 
was first discovered by |Carrington ( |1858[ ) and is beautifully shown in the Maunder 



butterfly diagram in [Figure 1 That sunspots harbor strong magnetic fields led to the 
discovery that the polarity of sunspot pairs is the same in a given solar hemisphere and 
opposite in the other during the cycle, and that it reverses from one cycle to the next 



( Hale et aL||1919[ ). Despite more than a century and a half of study, the solar dynamo 



remains poorly understood (for a review see |Hathaway|2010| ). 
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Proxy-magnetometry has also contributed substantially to the understanding of 
small-scale magnetic fields. 'Hale (1922 1 noticed that if spots do not appear in bipolar 
form, they are trailed or preceded by regions permeated by small, localized brightenings 
called 'faculae'. He found that magnetic fields are also present in those regions, and 
theorized about the existence of 'invisible sunspots'. [Sheeley ( |1967 1 first observed 
small magnetic elements well away from sunspots in the Can network. [Beckers &^ 



Schrdter ( 1968 1 observed the same features in spectra of plage and noted that they were 



abundant in the vicinity of sunspots, but much more rare in quiet sun. Observations 
by Mehltretter^ ( .1974) had adequate resolution to resolve the network into strings of 
small "bright points" located in intergranular lanes, that change in shape and size on 
timescales comparable to the lifetime of granules. The network had previously been 
associated with kilogauss field by (Stenflo,1973^ , suggesting that many of the observed 
structures were related. Direct evidence that faculae, gaps, and bright points were all 
manifestations of the same phenomenon was eventually provided by high-resolution 
observations, simultaneous in multiple wavelengths, of both disk and limb targets in 
a study by [Wilson ( 1981 1. The bright points form a dense pattern in plage and active 
network, while outside of active regions, they clump in patches that partially outline 
supergranular cells. 

The relative ease of observing magnetic fields via proxy-magnetometry resulted in 
extensive studies of bright points in the photospheric network. Imaging in wide-band 
CaiiH and K or Hor were the diagnostics of choice, until Muller & Roudier (1984 1 
switched to using the Fraunhofer G band around 430.8 nm in order to reduce the effects 
of chromatism in their telescope. The blue wavelength of the G band allows for both 
high resolution when compared to H o- as well as good sensitivity of film and silicon 
detectors when compared to Ca ii H and K. It is still widely used as one of the principal 
diagnostics for proxy-magnetometry (e.g., Rouppe van der Voort et al. 2005 1 Chitta 
et al.|2012t and references in Sect. 5 of 'De Wijn e t al.|2009| ). 



The suggestion by [Spruit & Zwaan ( 1981 ) that facular brightness enhancement 
is the result of radiation escape from hot granular walls was confirmed by 3D MHD 
simulation by Keller et al. (2004 1 and Carlsson et al. ( 2004| ). They modeled flux tubes, 
then "observed" them as if close to the limb using radiative-transfer codes to calculate 
the emergent intensity. Partial evacuation of the flux tube allows the observer to look 
deeper into the hot granular wall than would be possible if the flux tube were absent. 



2.2. Magnetograms 

The next step up from proxy-magnetometry is the "magnetogram", occupying the space 
between a qualitative diagnostic like G-band imagery and a quantitative diagnostic such 
as spectro-polarimetric inversions. Magnetograms are made using a narrow-band filter 
combined with a simple polarimeter. By tuning the filter to the wing of a magnetically 
sensitive line and observing circular polarization, a measure of the line-of-sight mag- 
netic field is found. While the measurement is not truly quantitative, it can be calibrated 
with some accuracy to provide magnetic flux density, i.e., the product of the filling 
factor and strength of the magnetic field. Because magnetographs are only sensitive 
to circular polarization, they only provide a diagnostic to the line-of-sight component 



of the magnetic field. [Figure 2| displays an example of a magnetogram together with 
proxy-magnetometry using imaging in the G band and in the Ca n H line. 

Perhaps the best known example of a magnetograph is the Michelson Doppler Im- 



ager (MDI, Scherrer et al.||1993j ) on the Solar and Heliospheric Observatory (SOHO, 
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Figure 2. Sample network region. Left: calibrated Fei 630.2 nm line-of-sight 
magnetogram, scaled between -1 X 10^ (black) and 10^ Mxcm"^ (white). Middle: 
G-band intensity. Right: Ca ii H intensity. Bright points in the G band and in the 
Ca n H line images correlate well with positions of concentrated field. The network 
consists of strings of several adjacent bright points, located in the intergranular lanes. 
Bright points in the Ca ii H image appear more extended and fuzzy than in the G- 
band image due to expansion of the flux tube with height. These images were taken 
with Hinode/SOT on March 30, 2007, ai'ound 00:24:3 UT. The coordinates indicate 
distance from sun center, so that fj. ^ 0.88. From ,De Wijn et al. ( |2009) . 



Domingo et al.|1995 1. The uninterrupted, long-duration, and seeing-free MDI observa- 
tions led to a huge increase in the understanding of, e.g., the formation of the photo- 
spheric magnetic network ( Schrijver et al.|1997 1 and dispersal of magnetic flux (Hage- 
naaretal.|1999l ). 



2.3. Quantitative Measurement 

The advent of the CCD camera for the first time allowed direct digitization of observa- 
tions and consequent analysis by computers, and as computers became more powerful it 
became feasible to compute spectral lines based on a polarized radiative transfer calcu- 
lation. This development led to the process now commonly referred to as "inversions". 
It is described in detail in 



section 3 



Simply put, from an observation of a spectral line 
one can derive the plasma parameters in the region of the solar atmosphere where that 
hne formed. 

The instrument of choice for these kinds of measurements is known as a spectro- 
polarimeter: the combination of a spectrograph and a full-Stokes polarimeter. The 
typical layout of such an instrument consists of the telescope feed optics, followed by 
a polarimetric modulator and a tradition spectrograph. A polarimetric modulator is a 
device that alters the polarized light entering the telescope in a known way. Because 
detectors are sensitive to intensity only, an 'analyzer', typically a linear polarizer, is 
then used to imprint polarization information into intensity. Out of at least four inde- 
pendent measurements one can reconstruct the Stokes vector consisting of intensity (I), 
linear polarization (Q and U), and circular polarization (V). An example of a spectro- 
polarimetric observation, taken with the Hinode/SOT Spectro-Polarimeter, is shown in 



Figure 3 The spectro-polarimeter of Hale ( 1908 1 could only modulate Stokes V, and 
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as a result it was only capable of line-of-sight diagnostics. While modern spectro- 
polarimeters are significantly more sophisticated than that instrument, the principle of 
operation remains the same. 

The first data to ever be processed with an inversion technique was taken with the 
Stokes I spectro-polarimeter developed at the High Altitude Observatory (HAO) and 
deployed at the Dunn Solar Telescope at Sac Peak. The success of Stokes I and its 
successor Stokes II led to the development of the Advanced Stokes Polarimeter (ASP, 
Elmore et al.||1992 i and many other instruments in the early 1990s, giving scientists a 



new window on magnetism in the solar photosphere (e.g., |Lites et al.|1987 1993j). As 
an example, the SOT Spectro-Polarimeter on Hinode is a derivative of ASP. 

Narrow-band imaging instruments have come a long way since the first birefrin- 
gent filters were built by |Lyot| ( [r944.) . Some instruments have sufficient spectral reso- 
lution that they can be used as true spectro-polarimeters. Fabry-Perot interferometers 
([Fabry & Perot ' 1 902) in particular have been very successful (e.g., |Scharmer et al.|2008j 



Viticchie et al. 2009p|Beck et al.|2010t Judge et al. 2010). The same computational tech- 



niques used with data from spectrographs for the inversion of plasma parameters can be 
applied to data from these instruments. The benefit of an imaging spectro-polarimeter 
is that it samples a 2D region, and usually takes only a few seconds to make a spectral 
scan. A slit-scanning spectro-polarimeter must raster to build up an image, a process 
that can take many minutes. The slow scanning is traded against high spectral resolu- 
tion and high sensitivity. 



3. Inversions 



Figure 4 shows a schematic overview of the process of spectro-polarimetric inversions. 
While this process is called an inversion, it really involves solving the forward problem 
of synthesizing an observation based on an atmospheric model and external fields, and 
comparing that synthesis to the actual observations. The fields and the model are then 
adjusted until a match is found. 

An atmospheric model is chosen and a guess is made of the external parameters. 
The radiation output can be calculated through radiative transfer if the atomic excitation 
is known. The atomic excitation, however, relies on a statistical equilibrium calcula- 
tion that depends on the radiation field. This loop must be iterated until convergence is 
reached. Then, the radiation output can be compared with the observations. If the qual- 
ity of the match is found to be insufficient, a new set of parameters must be evaluated, 
until a set is found that provides a good match. 

Many photospheric lines are formed in local thermal equilibrium (LTE), and statis- 
tical equilibrium does not need to be considered. Atomic excitation can be calculated 
directly from the plasma parameters. The calculation of the radiation output is fast 
because the consistency loop is broken. The simplest so-called 'Milne-Eddington' at- 
mosphere, in which the continuum and line source functions are the same and radiative 



transfer is easy to compute, is commonly used. Not surprisingly a method by Auer 



et al. ( 1977 1 that used this approximation was the basis for the first applications of 



inversions of spectro-polarimetric observations ( Lites & Skumanich|1982 ; Skumanich 



& Lites 1987V It remains the workhorse for inversions of photospheric vector mag 



netic field (Borrero et al.|20lT] ). More sophisticated LTE atmospheres may, e.g., derive 
height dependent information. They have more parameters than the Milne-Eddington 
atmosphere and are consequently computationally more expensive. 
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Figure 3. Example of spectro-polarimetric data taken by the Hinode/SOT 
Spectro-Polarimeter instrument on Dec 12 2006 at 11:50:04 UT. The X position 
of the slit was 242.5" W of disk center. From left to right: intensity (Stokes I), 
linear polarization (Stokes Q and U), and circular polarization (Stokes V). Around 
Y = -75" the spectrograph slit crosses a sunspot. The Zeeman effect causes the 
intensity profiles to be broadened in the Fe i 630.15 nm line and separated into three 
components in the Fe i 630.25 nm. The blue and red components are circularly po- 
larized in opposite directions, as can be seen in the Stokes V image. 



For lines formed in non-LTE (e.g., those that sample the chromosphere), these 
approaches are not sufficient. The full process has to be evaluated, requiring knowl- 
edge of the radiation field as well as a time-consuming iteration of the loop to ensure 
consistency. However, in some cases it can be assumed that the single-scattering ap- 
proximation is valid. If so, the statistical equilibrium calculation does not require the 
radiation output, breaking the loop and making the problem tractable. 

The codes that perform these calculations must find a global minimum in multi- 
dimensional parameter space that usually has many local minima. For the simpler at- 
mospheric models, i.e., Milne -Eddington and those using the assumption of LTE, the 
code will usually start the optimization using a Monte-Carlo method or genetic opti- 
mization in order to find the global minimum. As the solution improves, the code will 
switch to a derivative-based method to finalize the optimization. An example of an 



inversion is shown in Figure 5 



For the non-LTE case, even under the assumption of single-scattering, the calcula- 
tions are typically too costly to use a brute-force approach to find the global minimum, 
and a smarter approach is called for. Codes have been developed that find the solu- 
tion through a technique known as principle component analysis. A large database of 
lines is computed for many different parameter sets. This database is then decomposed 
into a basis of components called principle components. Once the database has been 
computed and the basis has been established, an observed line can be quickly approx- 
imated by a linear combination of elements of this basis, thus yielding an inversion of 
the observation. This method has been applied successfully to spectro-polarimetric ob- 
servations of chromospheric structures such as solar prominences ( Casini et al.|2005 1, 
and active region filaments ( Kuckein et al.|20()9 l. 
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Figure 4. Schematic overview of the process of spectro-polarimetric data inver- 
sion. In dark grey are inputs to the process: the observations that are to be inverted, 
the atmospheric model chosen by the scientist, and the external fields. In white are 
outputs generated by the radiative transfer and statistical equilibrium calculations 
(light gi-ey). 



Regardless of the inversion strategy, the nature of the data leaves us with ambigui- 
ties. The 180° ambiguity in the determination of the azimuth is well-known. A number 
of manual and automatic processes have been developed that attempt to remove this 
ambiguity by making assumptions about the magnetic field configuration that typically 
involves minimizing some quantity that depends on the azimuth of the magnetic field, 
e.g., the vertical gradient of the magnetic pressure. However, it is known that at least in 
some areas of the Sun these quantities are not in a minimum state, and so the resolution 
cannot be expected to yield the correct result. A good overview of this problem and the 
performance of several algorithms is given in Metcalf et al. (2006 1. 



4. Future Developments 

Many instruments now regularly provide us with a measure of the vector magnetic field 
in the solar photosphere. The vector magnetic field in the photosphere is determined 
a few times per hour from full-disk SDO/HMI data. While such measurements have 
greatly improved our understanding of the solar atmosphere and Sun-Earth interactions, 
it is clear that they are not sufficient to resolve a number of long-standing problems, 
such as coronal heating, the origin of the solar wind, and acceleration of coronal mass 
ejections. The key to all these unresolved issues lies with the magnetic field in the 
chromosphere and the corona. 

One can attempt to extrapolate the photospheric field measurements into the chro- 
mosphere and corona. Unfortunately, this process is difficult (see, e.g.,[De Rosa et al. 
2009[ ). Potential field extrapolations are clearly incorrect. More sophisticated methods 
have been applied with some success, but rely on "preprocessing" of the photospheric 
magnetic field in such a way that the extrapolated field approximates the morphological 
structure observed in a diagnostic of the chromosphere (typically Ho' or IR Can). A 
direct measure of the chromospheric vector magnetic and velocity field over the full 
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Figure 5. Example of an inversion of spectro-polarimetric data taken by the Hin- 
ode/SOT Spectro-Polarimeter i nstrument on Dec 12 2006 from 10:10:08 UT to 
13:11:18 UT. The data shown in Figure 3 was taken from this map. Clockwise from 
the top left (black to white scaling): continuum intensity (arbitrary), line-of-sight 
component of the magnetic field (-2 to 2kG), transverse component (0 to 2kG), 
azimuth (0° to 180°, i.e., with 180° ambiguity), incHnation angle (0° to 180°), and 
field strength (0 to 3 kG). 



solar disk is highly desirable to forgo this step, as it provides a true force-free boundary 
at the base of the heliosphere. 

The first successful measurements of coronal magnetic field have been made by 



Lin et al. ( 2004 1 using a fiber-fed integral-field spectro-polarimeter and by Tomczyk 



et al. (2008) using the Coronal Multi-channel Polarimeter (CoMP). All coronal field 



measurements by necessity employ a coronagraph that blocks out the light from the 
solar disk. Consequently it is only possible to determine the magnetic field integrated 
with some weighing along the line of sight in the plane of the sky. 

A number of instruments have been or are currently being developed with the 
intent of measuring magnetic field in the chromosphere, focusing on four lines that 
are thought to be the most promising: the IR Can triplet and Hei 1083.0 nm. The 
latter is easier to interpret and has been used with some success to measure magnetic 
field in prominences, active region filaments, and spicules ( Casini et al.|2005; Kuckein] 
et al. 2009; Centeno et al.||2010 ). However, in quiet sun the hne can be three orders 



of magnitude weaker than in active regions. The IR Can triplet at 849.8 nm, 854.2 nm 
and 866.2 nm is strong everywhere, but suffers in other aspects: the radiative transfer 
calculations required to model it are difficult, it forms over a broad height range, the 
lines are spectrally far apart, and 849.8 nm and 866.2 nm have blends. However, these 
four lines are still the most promising ones in the visible and near-infrared spectrum. 

The Coronal Solar Magnetism Observatory (COSMO) has been proposed as a fa- 
cility for synoptic measurements of the chromospheric and coronal magnetic field. It 
is to contain a suite of three instruments: a white-light coronagraph for coronal density 
measurements (under construction), an instrument for chromospheric magnetic field di- 
agnostics (under development), and a large coronagraph for coronal field measurements 
(in preliminary design phase). 
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It is also evident from theoretical considerations and simulation alike (e.g., |Rempel 



et al. 2009) 1 that high-resolution magnetometry is another frontier that must be explored 



Here, we have reached the limit of present-day telescopes. In the near future results are 
expected from the 1.6 m New Solar Telescope at the Big Bear Solar Observatory and the 
1.5 m GREGOR telescope on Tenerife, Spain. Under construction and development, 
respectively, are the Advanced Technology Solar Telescope and the European Solar 
Telescope. Both will have 4 m apertures, are designed with polarimetry in mind, and 
will be outfit with imaging and slit-scanning spectro-polarimeters. 

Acknowledgments. The National Center for Atmospheric Research is sponsored 
by the National Science Foundation. 
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